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Abstract 



We investigate the feasibility of probing the neutrino mass hierarchy and 
the mixing angle 1^13 with the neutrino burst from a future supernova. An 
inverse power-law density p ~ r n with varying n is adopted in the analysis as 
the density profile of a typical core-collapse supernova. The survival probabil- 
ities of v e and v e are shown to reduce to two-dimensional functions of n and 
013. It is found that in the n — sin 2 </>i3 parameter space, the 3D plots of the 
probability functions exhibit highly non-trivial structures that are sensitive to 
the mass hierarchy, the mixing angle c/>i3, and the value of n. The conditions 
that lead to observable differences in the 3D plots are established. With the 
uncertainty of n considered, a qualitative analysis of the Earth matter effect 
is also included. 
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I. INTRODUCTION 



A better knowledge of the neutrino property is regarded as one of the crucial keys in 
searching for new physics beyond the standard model. Recent experiments of neutrino 
oscillation have been able to uncover part of the neutrino properties, such as the neutrino 
mixing angles and the mass squared differences responsible for the solar neutrino problem 
and the atmospheric neutrino deficit [1-5]. However, there exists an ambiguity in the sign of 
the squared difference Am 32 = m§ — m 2 involved in the oscillation of atmospheric neutrinos. 
The two scenarios are referred to as the normal (m 3 ^> m 2 , mi) and the inverted (m 2 , m-i 
m 3 ) mass hierarchies. Furthermore, although an upper bound on the mixing angle 13 is 
established by sin 2 i3 < 0.02 [6,7], a definite determination of this angle has not been 
achieved yet. 

The rich physical content of a core-collapse supernova makes the supernova neutrino one 
of the most promising tools for the study of unknown neutrino properties and the supernova 
mechanism [8]. The supernova neutrinos are unique in that both neutrinos and antineutri- 
nos are produced at very high densities and high temperatures before propagating through 
matter of varying densities. Due to the wide range of matter density in a supernova, the 
neutrinos may go through two or even three (if the regeneration effect due to the Earth 
matter occurs) separate flavor conversions before reaching the terrestrial detectors. Further- 
more, the matter-enhanced oscillations [9,10] in a core-collapse supernova lead to a striking 
feature that a small variation of the mixing angle 13 can significantly alter the neutrino 
spectra. For supernova neutrinos, the main physical consequence arising from the ambiguity 
of the mass hierarchy is that both the higher and the lower level crossings occur in the v 
sector if the mass hierarchy is normal, while the higher crossing occurs in the v sector and 
the lower crossing occurs in the v sector if the mass hierarchy is inverted. 

The future galactic supernova is capable of inducing roughly 10 4 neutrino events at 
the terrestrial detectors, and is expected to provide a much better statistics than the 
SN1987A [11] did. This promising characteristic has motivated a wealth of discussions on 
how the neutrino fluxes from a supernova can facilitate the search of the unknown neutrino 
properties [12-18]. As generally realized, the main difficulty in extracting information from 
the supernova neutrinos arises from the poorly known exploding mechanism. Incomplete 
knowledge of the supernova leads to, among others, an uncertainty in the density profile of 
a supernova. 

The supernova neutrinos are usually assumed to propagate outward through an inverse 
power-law matter density, p ~ r n , with n = —3. However, due to lack of statistically 
significant real data, there is no clear evidence showing that the density distribution p ~ r -3 
provides a reliable connection between the dynamics of flavor conversion and the expected 
neutrino events at the detectors. In addition, the shock propagation in a supernova [19] 
represents a time-dependent disturbance to the matter density and causes a sizable effect 
to the neutrino flavor conversion. Since only the matter density near a resonance point is 
relevant to the flavor conversion and any local deviation from n = —3 cannot be ruled out, 
the profile p ~ r~ 3 should not be considered as a satisfactory description to the density 
shape for the purpose of extracting neutrino properties from the observation of supernova 
neutrinos. 

In this work, possible consequences resulting from variation of the density profile are 
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examined. The aim is to analyze the neutrino survival probabilities and to study how 
the uncertainty in n would affect the determination of the mixing angle 0i3 and the mass 
hierarchy. This paper is organized as follows. In Section II, we outline the n-dependent 
formulation of the survival probabilities for v e and u e . In section III, parameters obtained 
from the solar, the atmospheric, and the terrestrial experiments are taken as the input for 
constructing the 3D plots of the probability functions. The probability functions under 
both the normal and the inverted mass hierarchies are analyzed. In section IV, we discuss 
the feasibility of probing the neutrino mass hierarchy and 13 with the properties of the 
probability functions. In section V, the discussion is expanded to include the Earth matter 
effect. We then summarize this work in section VI. 



II. CONVERSION PROBABILITIES 

The density profile of matter encountered by the propagating neutrinos plays a crucial 
role in the dynamics of flavor conversion. In the literature, the neutrino flavor conversion 
in media of various density profiles has been widely discussed. However, the exact solution 
is obtained only for a few specific density distributions : the linear, exponential, hyperbolic 
tangent, and the 1/r profiles. It was suggested [20] that for an arbitrary inverse power-law 
density p ~ r n , an extra correction factor F (a function of n and the mixing angle) can be 
supplemented to the standard Landau-Zener [21] formulation of the level crossing to account 
for the effect due to deviation from a linear density profile. 

With the extremely high electron number density in a supernova, the effective mixing 
angles in matter for the neutrino and the antineutrino become 0^ ~ n/2, 9^ pa n/2 and 
0™ pa 0, ^2 ~ 0) respectively. Using the standard parametrization of the neutrino mixing 
matrix: U 2 X = cos 2 0i 3 cos 2 9i 2 , U 2 2 = cos 2 0i 3 sin 2 6\ 2) and U 2 3 = sin 2 0i 3 , the survival 
probabilities for v e and v e can be written respectively as [22,13] 

Prior = U 2 el P t P h + [/* (1 - P t )P h + C/ 2 3 (l - P h ), (1) 

Pnor = U 2 el (l-P l ) + U 2 e2 P l , (2) 

for the normal hierarchy and 

P tnv = U 2 e2 (l-Pi) + U 2 el P h (3) 



Pinv = U 2 e2 P x P h + U 2 el (l - P l )P h + C/ 2 3 (l - P h ), (4) 

for the inverted hierarchy. Note that Ph {Ph) and Pi (Pi) represent the higher and the lower 
level crossing probabilities for u e (z/ e ), respectively. For arbitrary density profile and mixing 
angle, the Landau-Zener formula is modified as 

exp(-f7^j - exp(-f 7MjJ%) 
1 - exp(-^h,i—TQ-) 
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where 7^ are the adiabaticity parameters, 0^ are the mixing angles between the ith and 
the jth mass eigenstates, and F^i are the correction factors to a non-linear profile. Note 
that Ph — Phi and that P\ can be obtained directly from P\ by replacing 0jj with tt/2 — G^-. 

For a typical core-collapse supernova, the electron number density can be written as 
N e = (^-)cr n , where Y e is the electron number per baryon, m n is the baryon mass, and 
c is a constant representing the scale of the density 1 . The adiabaticity parameter for this 
density profile has the general form 

lKl 2|n| l E } l cos2e,/ l 2 v / 2G F ^ C J ' K) 

where Gf is the Fermi constant, E is the neutrino energy, and Am?- = |m 2 — m 2 |. In the 
numerical calculation, it would be more convenient to write Fh and F as the Euler integral 
representation of the hypergeometric function: 



2n- 1 



2n )^( 2 2n ) 

The above expressions for P h j, 7^, and F/^ can be applied to an arbitrary profile and to 
both large or small mixing angles. It was pointed out [24,25] that there exists a subtlety in 
the physical meaning of resonance conversion: For the large mixing angle, the adiabaticity 
parameters, 7^ and 7?, should each be calculated at the point of maximum violation of 
adiabaticity (PMVA) instead of the point of resonance. Note that while the values of 7^ 
depend on the locations where they are calculated, the values of jhFh an d J1F1 remain 
invariant. To simplify the calculation, we choose to evaluate jh^h and 7/i 7 } at the locations 
of resonance. 



III. NUMERICAL ANALYSIS OF P AND P 

With the numerical input of Am^i, Am^, 6>i 2 , and c (see, e.g., ref. [14] and the refer- 
ences therein), the survival probability functions P = P(Am 2 1 , Am^, 13 , 9 12 , E u , n, c) and 
P = P^Am^, Aml 2 , (fi 13 , 9 12, E p ,n,c) reduce to P = P(E U , n, 0i 3 ) and P = P(E p ,n,(f)is), 
respectively. The ambiguity of the mass hierarchy gives rise to four distinct probability func- 
tions to be investigated: P nor , Pnor for the normal hierarchy and Pi nv , Pi nv for the inverted 
hierarchy. Their properties can be examined by a series of 3D plots in the E — sin 2 i3 
parameter space. As an illustration, we show only the 3D plots of P nor for several n values 
in Fig. 1. From the 3D plots of P nor , P nor , Pi nv , and Pi nv , it can be concluded that all 
the probabilities exhibit no significant energy dependence except for the low energy end. 
This behavior implies that when the adiabaticity parameters in Eq. (6) are calculated, the 
impact coming from the variation of n and 0i 3 dominate over that of the variation of energy 
in the typical range, E < 10 2 MeV. 



The value of c varies very weakly with r over the range 10 12 g/cm 3 < p < 10 5 g/cm 3 [23]. 
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Since the neutrino population is extremely small at the low energy end of the spectrum, it 
would be convenient to simply take the average energies, e.g., (E u ) = 12 MeV and [E v ) = 15 
MeV, in the numerical calculation. This approximation further reduces the probabilities to 
functions of only n and 0i 3 . The survival probabilities for u e and u e can be plotted in the 
n — sin 2 0i3 space, as shown in Fig. 2 and Fig. 3, respectively. It is seen that if n < —6, 
the values of all the probability functions approach a constant ~ 0.6 regardless of the mass 
hierarchy and the value of 0i3. Thus, the mass hierarchies are indistinguishable and the 
information about i3 is lost if n < —6. As n increases from —6, the probability functions 
of u e in Fig. 2 are seen to drop through a transition near n ~ —5, while that of u e in Fig. 
3 jump through a transition near n ~ —4.5. Furthermore, P nor and P inv exhibit an extra 
non-trivial structure for n > —4. In the following discussion, we divide n into three regions: 
n < —6, — 6 < n < —4, and n > —4. The probability functions for u e and u e shall be 
discussed separately. 



A. -Pnor Pinv 



We first note that the condition P nor = P inv is satisfied if the higher crossing is ex- 
tremely non-adiabatic: Ph — > 1, as implied by Eqs. (1) and (3). Thus, P nor and Pi nv are 
indistinguishable if the values of n and sin 2 13 result in Ph — > 1, which occurs in part of the 
n — sin 2 0i3 parameter space, as can be seen by comparing Figs. 2(a) and 2(b). To account 
for the n-dependent transition of both P nor and P inv in the range — 6 < n < —4, we note 
that the adiabatic parameter at the lower crossing takes the form 

0.43 x 10" 5 , o nx i 
7/ ^ r-j (0.39 x 10~ 30 )", (8) 

which yields 7j C 1 and Pi ~ cos 2 6i 2 (non-adiabatic transition) for n < —6. On the 
other hand, the same expression leads to 7z ^> 1 and Pi m (adiabatic transition) when 
n > —4. It is clear that 7; goes through a transition from 7; 1 to 7^ ^> 1 as n varies from 
n ~ — 6 to n ~ —4. A simple calculation shows that P inv ~ P nor « 0.6 at n = —6, while 
P inv ~ P nor fa 0.3 at n = — 4 if sin 2 i3 < 10~ 3 . This result implies that the uncertainty of 
n between n = —4 and n = —6 could lead to a variation of the survival probability by a 
factor of two and complicate the interpretation of the v e events. 

The 3D plots of Pi nv and P nor in the n — sin 2 i3 space become distinguishable for n > — 4 
if Ph 7^ 1. Note that for n > —4, Pi fa and P nor becomes 

P nor ~ sin 2 0i3 + Ph(sin 2 6 l2 cos 2 013 - sin 2 i3 ), (9) 

where Ph is given by Eq. (5), and the adiabaticity parameter 7^ in Ph is given by Eq. (6): 

lh * ^(37.6 x 10- 30 x cos2#^^H. (10) 
\n\ cos 20i3 

For n > —4, the subtle dependence of P nor on n and sin 2 i3 can be seen clearly from 
Eq. (9) and Fig. 2. For small \n\ and large 0i3 (near the right corner of Fig. 2(a)), the 
higher level crossing is adiabatic (7^ ^> 1 and P h ~ 0). Thus, the first term in Eq. (9) 
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dominates and P nor ~ sin 2 0i 3 <C 1. On the other hand, the higher level crossing becomes 
non-adiabatic (7^ <C 1 and Ph ~ 1) for relatively larger |n| and smaller 13 . Thus, the 
second term in Eq. (9) begins to dominate, and P nor ~ sin 2 9 12 ~ 0.3. Note that since the 
two values of P h , P h « 1 and P ft « 0, give rise to the above two distinct values of P nor 
representing the two sides of the fast transition area in Fig. 2(a), the condition P h — 1/2 
should reasonably describe the fast transition of P nor for n > —4. Furthermore, due to the 
smallness of the upper bound on 0i 3 , the arguments in the numerator of Eq. (5) satisfy the 
relation 

71 TP ^ 71 ^hFh ni x 
-ZlhFh < TT • 2 7 , 11 

2 2 sin 0i3 

which implies that 

exp(-| 7fc F h ) » 1 (12) 
at the narrow transition region. The condition Eq. (12) leads to G h m 1, where 
r 71 1 ( Am 32 w llz CQS 2 ^3 1 sin 2 20 13 

Gft = 4(hT2) R ( ^T )( c ) " ( ^2^ ) ^ (13) 

Since ~ 1 in the region of interest (n > —4 and sin 2 0i 3 < 1CT 2 ), the above condition 
can be approximated as 

G h {nAn) * ^ ^°" 4 (37.6 x 1(T 30 x cos20 13 )^^^ » 1. (14) 
\n\ cos20i 3 

It can be shown that Ph ~ 1 if Gh{n, 0i 3 ) < 1 and if G/ l (n, 0i 3 ) > 1. Take n = —3 for 
example, the sudden probability transition is located near sin 2 13 ~ 10~ 5 . Thus, P ~ if 
sin 2 13 > 10 -5 , which is unique to the normal mass hierarchy. However, if sin 2 13 < 10~ 5 , 
then P ~ 0.3 for both the normal and the inverted mass hierarchies. This feature is clearly 
seen in Fig. 2. 

Eq. (14) suggests that a slight variation of the power may cause an ambiguity in the 
interpretation of 0i 3 and the mass hierarchy that are derived from the observation of neutrino 
events. Note that although the numerical values in Eq. (14) vary with the input parameters, 
the physical content remains unaltered. We summarize this subsection as follows. 

1. Given the input values of Am 2 1? sin 2 # 12 , and c, it can be shown that Pi is adiabatic 
(7i ^> 1) for n > — 4 and non-adiabatic (7/ 1) for n < — 6. 

2. For the normal mass hierarchy, the two distinct values of P due to Ph ~ 1 and Ph ~ 
for n > —4 are separated by the condition Eq. (14). 

3. In principle, a direct observation of the v e events could be used to distinguish the mass 
hierarchy if the values of n and i3 result in P h ^ 1, as suggested by Figs. 2(a) and 
2(b). 
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E. Prior P%nv 



The survival probabilities of v e are described by Eq. (2) and Eq. (4) for the normal and 
the inverted mass hierarchies, respectively. The 3D plots of P nor and P inv are shown in Fig. 
3. Note that P h = P h) and that 7/, F h and P t can be obtained respectively from 7^, F h and 
Pi by the swap sin0i 2 <-> cos0i 2 . Since P; ~ sin 2 9 12 ~ 0.3 and P^ ps 1 for n < —6, it follows 
that 

Pnor = Pinv — cos 2 13 cos 2 i2 (l - sin 2 i2 ) + cos 2 i3 sin 4 6>i2 w 0.6. (15) 

This explains why the mass hierarchies are also indistinguishable from observing the v e 
events if n < —6. In addition, the transition behavior of P nor and Pi nv for — 6 < n < —4 can 
be explained in the way similar to that of P nor and Pi nv . 

For n > —4, the lower level crossing becomes adiabatic: P t ps 0. If the higher crossing 
remains non-adiabatic (i\ Pi 1), it leads to P Ttor ~ Pj ni) ~ cos 2 0i 3 cos 2 0i 2 ~ 0.7, which is 
slightly higher than that for n < —6. However, when Ph departs from unity, there would be a 
sudden drop of the probability function if the mass hierarchy is inverted: Pi nv ~ sin 2 i3 <C 1, 
as shown in Fig. 3(b). The sudden drop of this probability function is similar to that of 
P ri ar for the neutrino sector, and can be characterized by the same condition for that of the 
neutrino, Eq. (14), with a slight change of the numerical values. The properties of all the 
probability functions for n > —4 and n < — 6 are summarized in Table 1. 

IV. POWER-LAW DENSITY PROFILE, 13 , AND MASS HIERARCHIES 

We shall now investigate whether and how the information about the mass hierarchy 
and 0i 3 could be extracted from the observation of supernova neutrinos. Although the 
uncertainty in the density profile is unavoidable due to incomplete knowledge of the super- 
nova mechanism, the results in Fig. 2 and Fig. 3 may still provide a useful guideline. We 
summarize the hints as follows. 

1. As suggested by Figs. 2(a) and 2(b), the mass hierarchy may be identified as normal 
if the survival probability of v e is observed to be extremely small, i.e., if P 0.3. On 
the other hand, an extremely small survival probability for 9 e , P 0.7, would signal 
an inverted hierarchy, as suggested by Figs. 3(a) and 3(b). Although this feature 
provides no information about the mixing angle 0i 3 , it predicts the mass hierarchy 
without knowing the details of the density shape, i.e., the exact value of n is irrelevant. 
Note that the numerical values would vary slightly with the input parameters. 

2. If P ~ P ^ 0.6 is observed, the power n of the density profile is limited to n < — 6 and 
a significant deviation from n = — 3 is implied. The information about mass hierarchy 
or is unavailable in this case. 

3. If P < P can be deduced from experiments, then n > —4, G h (n, 13 ) > 1, and the 
inverted mass hierarchy are implied. These conditions lead to sin 2 0i 3 > 4 x 10™ 4 . On 
the other hand, if P > P is observed, it implies that (i) the mass hierarchy is normal 
and 0i 3 is undetermined (Figs. 2(a) and 3(a)), or (ii) the mass hierarchy is inverted 
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and Gh{n,4>\z) < 1 (Figs. 2(b) and 3(b)). A further observation of the Earth matter 
effect may be useful in selecting the correct scenario. 

4. From Figs. 2(a) and 3(b), it is seen that 0i 3 can be sensitive to the direct observation 
of P nor or Pi nv in only part of the parameter space. As discussed earlier, the prediction 
of 13 depends crucially on the exact value of n in this part of parameter space. Thus, 
it is very difficult to establish a satisfactory constraint on 0i 3 from a direct observation 
of P or P alone. For example, Eq. (14) suggests that a deviation of ±0.5 from n = — 3 
could result in an uncertainty in the prediction of 0i 3 by up to two orders of magnitude. 
A better knowledge of n or an extended analysis that includes the Earth matter effect 
would help set the constraint of this tiny mixing angle. 



V. OBSERVATION OF EARTH MATTER EFFECTS 

The regeneration effect of the supernova neutrinos when crossing the Earth has been 
widely discussed [12-14,16,26,27]. Since the Earth matter density, qe ~ a few g/cm 3 , 
is roughly the same order of magnitude as the density for the lower level crossing in a 
supernova, the neutrino fluxes may receive sizable modification due to the oscillation effects 
in Earth. In general, the Earth matter effect is signaled by the flux difference observed at 
two terrestrial detectors and [12]: 

/(l) _ /e (2) „ p h{1 _ 2m pW _ p(2) )(/ _ ; ); (16) 



fP ~ /f * (1 - 2P;)(A ( e ) - P[?M ~ /°), (17) 

for the normal mass hierarchy, and 

/e (1) " fP ^ (1 - 2fi)(P£> - P 2 ( e 2) )(/ e ° - f° x ), (18) 



# - /f ~ P h (l - 2P{){PS - pg)(fS - /°), (19) 

for the inverted mass hierarchy. In the above expressions, (f^) is the observed v e (P e ) 
flux at the detector Pj l J {Pf^) is the probability that a Uj (pj) arriving at the Earth 
surface is detected as a u e (P e ) at the detector, and f® >x (/g 5 ) is the initial flux for the specific 
neutrino (antineutrino), with x = fi,T. 

The Earth matter effects could affect (i) v e flux only; (ii) v e flux only; (iii) both v e and v e 
fluxes. We assume that the suppression of matter effect, if any, is due solely to the smallness 
of Ph. Possible consequences of the above three scenarios are summarized in Table II and 
discussed below: 

1. If the Earth effect is observed only in the v e flux (/« - /j 2 ) ^ and ff ] - /j 2) = 0), 
it requires an extremely small Ph under the inverted hierarchy, as can be seen from 
Eqs. (18) and (19). Results in Table I show that Ph — > is possible when n > —4 and 
Gh{n, 0i 3 ) > 1. These conditions lead to a lower bound on 13 : sin 2 13 > 4 x 10~ 4 . 
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2. If the Earth effect is observed only in the u e flux (/« - = and tf ] - /J 2) ^ 0), 
then Eqs. (16) and (17) suggest that Ph is extremely small and the mass hierarchy is 
normal. This leads to the same constraint: sin 2 13 > 4 x 10~ 4 . 

3. If the Earth matter effect is observed in both the v e and the v e fluxes, then Ph ^ is 
required. It can be seen from Table I that this condition can be satisfied if (i) n > —4 
and Gh(n,(pi 3 ) < 1, or (ii) n < —6. Note that, as discussed earlier, the observation of 
supernova neutrino loses its predictive power if n < —6. Even from the first condition: 
n > —4 and Gh(n,(p 13 ) < 1, an evident constraint on 13 is still not available unless 
the uncertainty of n can be reduced significantly. Furthermore, this scenario provides 
no information about the mass hierarchy. 

There is a possibility that the constraint on 13 might be available from checking the 
signs of /j 1 -* — fP and /g — f'P , if the Earth matter effect is observed in both the v e and the 
v e fluxes. Suppose that one of the two detectors, D^ 2 \ is not shielded by the Earth matter. It 
then leads to the replacements: P^ — > \U e 2\ 2 and P^ — > |f/ e i| 2 , which satisfy > \U e 2\ 2 
and P$ > \U e i\ 2 . Since the average energies obey the hierarchy: (E(u e )} < {E(y x )) and 
(E(u e )) < (E(9 X )), there exists an energy E c (E c ) at which /? — /° (/? — /?) changes 
sign [12]. In general, /° - /° > (f° - f° > 0) in E < E c (E < E c ), and / e ° - /° < 
{fs ~~ fs < 0) m E > E c (E > E c ). Furthermore, the magnitude of Pi varies from (for 
n > —4) to cos 2 9 12 ~ 0.7 (for n < —6), and P\ varies from to sin 2 9\2 ~ 0.3. Thus, 1 — 2P\ 
is always positive, while 1 — 2Pi flips sign over the transition region — 6 < n < —4. The 
above arguments suggest the following: 

1. If — and /J 1 ^ — fg 2 ^ are both observed to be negative at the high energy end of 
the spectrum, or both positive at the low energy end, then it implies 1 — 2Pi > (from 
Table I), n > —4, and G h (n, 13 ) > I. This leads to the constraint: sin 2 13 > 4 x 10~ 4 . 
However, the mass hierarchy is undetermined from this result. 

2. If /W - fP and fP - /i 2) are of opposite signs, then 1 - 2P X < 0, and n < -6. No 
further information about 0i 3 or the mass hierarchy is available. 

We summarize the above results in Table III. 

VI. SUMMARY AND CONCLUSION 

The supernova neutrinos may provide a promising future for the study of unknown 
neutrino properties. However, the detailed knowledge of the core-collapse supernova event 
is still far from complete. In addition to the uncertainties in the original neutrino fluxes 
and in the effects due to the shock propagation, the original neutrino spectra can be further 
deformed by the flavor conversion when the neutrinos propagate through matter of uncertain 
density profiles. 

In this work, parameters obtained from recent experiments are taken as the input for 
the purpose of analyzing the survival probabilities of u e and V e . It is suggested that the 
influence coming from the energy variation can be excluded. The effort is then focused on 
investigating how the unknown mass hierarchy, the mixing angle 0i 3 , and the uncertainty 
in n would affect the probability functions. 
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It is shown that the non-trivial behavior of the probability functions can be well illus- 
trated by the 3D plots in the n — sin 2 13 parameter space, and that the uncertainty of n 
could lead to ambiguity in the interpretation of 0i 3 and the mass hierarchy. Roughly speak- 
ing, the probability functions behave differently in three regions of the parameter space: 
n < —6, — 6 < n < —4, and n > —4. As far as the mass hierarchy and the mixing angle 0i3 
are concerned, the information is lost if the supernova neutrinos encounter a relatively steep 
density profile (n < —6) near the location of flavor conversion. For a not as steep density 
profile (—6 < n < —4), all the probability functions go through a transition that is governed 
by the variation of n. This transition depends only very weakly on the mass hierarchy and 
0i 3 . For n > —4, the probability functions vary with the mass hierarchy, the value of n, 
and 0i3 in a non-trivial fashion, as depicted clearly by Figs. 2 and 3. Furthermore, this 
non-trivial structure is found to be divided by a function of n and 0i3 through the condition 
Eq. (14). 

For the qualitative observation of the Earth matter effect, it can be shown that the 
constraint on 13 would be available only if n > —4. However, the exact value of n is 
irrelevant to the constraint as long as n is greater than —4. 

It is hoped that Eq. (14) and the 3D plots of the probability functions could provide a 
guideline to finding useful observables from the future supernova neutrino experiments, and 
to better help shed light on the desired understanding of the neutrino properties. We shall 
return to this topic in the near future. 
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TABLES 





n > — 4 


n < 


—6 




Ph(Ph) 


0, if G ft (n, 0i 3 ) > 1 


2 

cos 


013 ~ 1 






1, if G h {n^i Z ) < 1 








Pi 





2 

cos 


6*12 ~ 0.7 




Pi 





. 2 

sin 


012 ~ 0.3 




P 


i "2/1 n i * 2 i 1 1 t~) \ 

normal: sin #12 -Hi + sin 0i3(l — Ph) 


■ 2 
sin 


f 12 + (COS 012 — Sin 012) cos 1 


1\2 ~ 0.6 




~ sin 2 0is k, 0, if Ghin, 6^) > 1 










~ sin 2 6»i2 « 0.3, if G h (n, i3 ) < 1 










inverted: sin 2 0u ~ 0.3 








P 


normal: cos 2 0i3 cos 2 6*12 ~ 0.7 


cos 2 


6>i2 + (sin 2 6>i2 - cos 2 12 ) sin 2 1 


9i2 ~ 0.6 




inverted: cos 2 0\2Ph + sin 2 4>is(l — Ph) 










~ sin 2 ^13 « 0, if G h (n, 0i 3 ) > 1 










~ cos 2 9i2 « 0.7, if Ghin, (f> 13 ) < 1 









TABLE I. Properties of the probability functions for n > — 4 and n < — 6. 



matter effect requirement prediction 



z/ e only 


Ph^O, 


n > -4, G h (n, 


013) > 1 


inverted, sin 2 i3 > 4 x 10 4 


v e only 


Ph^O, 


n > -4, G h (n, 


013) > 1 


normal, sin 2 13 > 4 x 10~ 4 


both f e and v e 




Ph^O 




* 



TABLE II. Predicting the mass hierarchy and ^13 from possible scenarios of the Earth matter 
effects. (* See Table III for further predictions.) 







/ 




requirement 


prediction 


E > E c , (E > E c ) 




+ 




/° < 0, 1 - 2Pi < 


n < — 6 


(high energy end) 






/e°- 


/° < 0, 1 - 2Pi > 


n > -4, sin 2 i3 > 4 x 10~ 4 


E < E c , (E < E c ) 


+ 


+ 


/e°- 


/° > 0, 1 - 2P { > 


n > -4, sin 2 i3 > 4 x 10~ 4 


(low energy end) 






fS- 


/° > 0, 1 - 2P t < 


n < -6 



TABLE III. Predicting 13 and n from the signs of / = / g (1) - /J 2) and / = fP - fP if the 



Earth matter effect is observed in both v e and u e fluxes. Note that the information about the mass 
hierarchy is unavailable. 
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FIGURES 






FIG. 1. The 3D plots of 

Pnor — Pnor{E u , (pis) for different values of n. The following values 
are adopted: Am§ 2 = 3.0 x 10~ 3 eV 2 , Am^ = 7.0 x 10~ 5 eV 2 , sin 2 12 = 0.8, and c = 7.0 x 10 31 g • 
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(a) normal hierarchy 




(b) inverted hierarchy 




FIG. 2. The 3D plots of P = P(n, <pis) under both (a) the normal, and (b) the inverted mass 
hierarchies. The average energy (E v ) = 12 MeV is adopted. 
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(a) normal hierarchy 




(b) inverted hierarchy 




FIG. 3. The 3D plots of P = P(n,4>is) under both (a) the normal and (b) the inverted mass 
hierarchies. The average energy (E p ) = 15 MeV is adopted. 
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